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Abstract. We apply the Separation of Variables method to obtain eigenvectors of commu- 
0^ ■ ting Hamiltonians in the quantum relativistic Toda chain at a root of unity with boundary 
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O ■ 1 Introduction 

During the last decade, a considerable progress in application of the Separation of Variables 
method to a broad class of quantum integrable models has been achieved. This progress was 
initiated by the paper [I| of Sklyanin who has proposed (using the i?-matrix formalism) a recipe 
for a Separation of Variables in the case of quantum Toda chain where the algebraic Bethe 
ansatz fails. The next important step was done by Kharchev and Lebedev [2] who realized an 
iterative procedure of obtaining the eigenfunctions of the n-particle open Toda chain by some 
integral transformation from the eigenfunctions of the (n — l)-particle open Toda chain. This 
iterative method was applied later to a relativistic Toda chain [3]. 

The Separation of Variables was found to be effective for the derivation of the eigenfunctions 
of the commuting Hamiltonians of the n-particle quantum Toda chain with boundary interac- 
tion [H [5] in the framework of Sklyanin's approach to the integrable boundary problems [6]. 
In the paper [5] the n-particle eigenfunctions of the quantum Toda chain when the first and 
last particles are exponentially interacting with the walls (the two-boundary interaction) is con- 
structed by means of an integral transformation of the eigenfunctions for the Toda chain with 
one-boundary interaction (the auxiliary problem). These eigenfunctions, in turn, are constructed 
using the eigenfunctions of the n-particle open Toda chain. 

Recently, a Separation of Variables was applied to the iV-state spin lattice model - the 
Baxter-Bazhanov-Stroganov (BBS) model [71 [8]. In the i?-matrix formulation of this model the 
cyclic L-operator [9j[TU] emerged. These operators are intertwined by the six- vertex i?-matrix at 
a root of unity. The cyclic L-operator of the BBS model reduces (at special values of parameters) 
to the L-operator of the quantum Relativistic Toda Chain (RTC) at a root of unity 
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The goal of this contribution is to find common eigenvectors of commuting Hamiltonians of 
the RTC at a root of unity with boundary interaction. We use for this the approach from [5]. 

Let us describe briefly the structure of the paper. We start with a short account of the model 
under consideration and remind its properties based on the presentation in [TT] (Section 2). We 
recall the standard Sklyanin formalism introducing a chain of L-operators satisfying some RLL = 
LLR relations as well as its "twisted" version (|2.3|) such that the corresponding ("twisted") 
monodromy matrix satisfies the "Reflection Equations" and some "unitarity" condition. We 
specify a "boundary matrix" solution K + (\) of the dual Reflection Equations which are necessary 
for a description of a two-boundary RTC. We close this section with a Hamiltonian generating 
function to the two-boundary RTC. 

We give the explicit formulas for the eigenvectors of the open RTC in the Section 3 using the 
eigenvectors formulas for the open Bazhanov-Stroganov quantum chain of which the open RTC 
is a (specified) particular case [9]. 

In Section 4 we consider the eigenvectors of the one-boundary RTC specifying some par- 
ticular choice of the boundary matrices. We identify the coefficients of the corresponding 
generating function with the commuting one-boundary RTC Hamiltonians and compute the 
common eigenvectors. The same questions are treated in Section 5 for the two-boundary RTC 
using the Separation of Variables method. We write the Baxter equation for separated variables 
explicitly (|5.3p . 

Section 6 contains an another proof of the commutativity of the one-boundary RTC Hamilto- 
nians of the Section 4. We show in fact that they belong to a wide class of Hamiltonians which 
can be expressed like a sort of Gelfand-Retakh quasideterminants. This quasideterminants are 
(in our case) just Pliicker coordinates A s Aq 1 on GV n _i(n). The commutation (and, hence the 
existence of the common spectrum) is guaranteed by the results of [12] (see also [131 114|). 

We finish with a short discussion in the Section 7 and we leave some tedious computations 
to the Appendix A. 

We are dedicating this paper to the memory of Vadim Kuznetsov. His important contributions 
to the subject of Separation of Variables and in particularly to its applications in quantum 
integrable models with boundary interaction [H fT5| fT6j [TT] are widely acknowledged in the 
Mathematical Physics community. In particular, it was Kuznetsov's paper |16] where the Se- 
paration of Variables was performed for the quantum relativistic Toda chain and commuting 
Hamiltonians of the quantum RTC with integrable boundary interaction (the quantum RTC for 
general classical root systems) were obtained. 

2 The formulation of the model 

In this section we recall briefly the subject of the model called the quantum Relativistic Toda 
Chain (RTC) chain at root of unity |llj . 

Let oj = e 2m / N , N > 2. For each particle k, k = 1,2, ... , n, of the n-particle RTC there is 
a corresponding iV-dimensional linear space (quantum space) V k with the basis | 7)^,7 G Zjy, 
and a pair of operators {u k ,v k } acting on V k by the formulas: 

Vk\l)k = vw 7 |7) fc , Mfc|7>fc = u\l ~ l)fc- 

The space of quantum states of RTC with n particles is V = Vi <8> V2 £S> ■ ■ ■ <8> V n . We extend 
the action of operators {u k ,v k } to V defining this action to be identical on V s , s 7^ k. Thus we 
have the following commutation relations 

U j U k = U k U jl VjV k = V k Vj, UjV k = LU 5 l' k V k Uj. 
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^ n (A) B n (X) 
C n {\) D n (X) 



Monodromy matrix for RTC with n particles is 
T n (X) = Li(A)L 2 (A) • • • L n (\) 
where 

L k {\) = 



1 + 



u k v k 







k = 1, 2, . . . , n. 



(2.1) 



Thus the model depends on three parameters u, v and x. It is possible to make this model to be 
inhomogeneous by attaching the index k to the parameters u, v and x. All the results presented 
in this paper can be directly generalized to inhomogeneous case. 

The straightforward calculation shows that the operators L k (X) satisfy the following commu- 
tation relations: 



3 the standard notations 
i?-matrix R(X, fj,) has the form 

R(\,fi) 



(2.2) 



where the standard notations (A) = L k (X) ® I, L? (//) = I (8) L k (fi) were used and the 



/ A - lju 

X-fi 

(l-w)A w(A-m) 

\ X-cufi J 

Using (|2.2p . we obtain the following relation for monodromy matrix 

R(X, M )T« (X)TP (jj) = ( M )r« (X)R(X, M ) . 

We will use the method [6] of Sklyanin for description of RTC with integrable boundary 
interaction. The L-operator (|2.1|) together with 



L fc (A) = Xa x L~ 1 (e/X)a x 



xX + oj9v k 1 

A«r 1 w fc 



Lj6>t> fc l a z L k {uj9 /X) T a z 



(2.3) 



satisfy the following relations 



fl(A,/i)L«(A)L^(/i) = L^(/,)L^(A)i?(A,/x), 
Z< 2) (n)R(X, p)L® (A) = 4 1} (X)R(X, n)L{ 2) 0*) , 
4 2) (fi)R(X, M )L« (A) = L « ( A)i?( A, M )4 2) (/,) , 



r(2) 



(2), 



where 



i2(A,/x) = naf> R{X,6 / ^)a, 



r (2) 

X 



<T 3 



1 

1 



f e-Xfi o 

w£ 



\ cjA/z — A/i 

1 
-1 



cj0 - \ 

A// 

lu6- Xfi 

w# — cjA/x / 
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Here 6 is an arbitrary parameter. It is important that R(fi, A) = R(X, fi) = P\2R(X, n)Pi2, where 
the operator P12 is the operator of permutation of two spaces. Then the matrix 

U(X) = L n (X) ■ ■ ■ L^K^L^X) ■ ■ ■ L n (\) = ( ^ *M ) 

satisfies the reflection equation 

R(X,n)U^(X)R(X,fi)U^(fi) = uW(n)R(\,v)U^(\)R(\,n) (2.4) 
if K_(X) satisfies the same equation: 

R(X, (X)R(X, ii)K^ (jj) = (lm)R(\, fi)K W (\)R(\, fi). (2.5) 

The solution K_(X) with non-operator matrix elements is 



a-(f-A) r?_+/?_f 
7?_+/3_A - A) 



This solution depends on four parameters a_, j3_, <5_, 77_. It can be found directly solving (|2.5p 
or adapting the formulas from [18] (see also [E]). The matrix £7(A) satisfies the so-called 
unitarity condition: 

"<^> -(«>«) 

1 / (A 2 - w 2 0)^4(A) A 2 (l-w)i3(A) + (^-A 2 )C(A) 



w (e _ A 2 ) V cu(u;# " A 2 )£(A) + w(9(l - w)C(A) (A 2 - w 2 0)P(A) 

We have ^4(A)^4(/i) = .A(/z).A(A) from the reflection equation f|2.4|> . It means that A(X) 
is generation function with respect to A for commuting set of operators. These operators are 
Hamiltonians for one-boundary RTC. 

To describe a two-boundary RTC we need the dual reflection equation: 

R(X, n)K® (ji)Q{\ fj) (A) = (A)Q(A, n)K® ( M ) J?(A, m) , 

where 

Q(A, M ) = ((^ 1 (A,M))" 1 ) tl . 
Its solution is 



K + (X) 



a+(^-A) r ?+ + /3 + f 

A 



??+ a; + /3 + A - u^A) 



Then the generation function for the Hamiltonians of the two-boundary RTC is given by t(A) 
Tt(K + (X)U(X)) or explicitly: 

t(A) = a + (9uj 2 /X - X)A(X) + {P+X + t] + uj)B{X) 

+ (P+9u/X + 77+)C(A) + 5 + {uj6/X - a>- 1 A)P(A). 
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3 Eigenvectors for the open RTC 

In order to give explicit formulas for the eigenvectors for RTC we remind the definition (see 
for example [19] ) of w p (s) which is an analogue of T^-function at root of unity. For any point 
p = (x, y) of Fermat curve x N + y N = 1, we define w p (s), s G Z^r, by 

WP(S) ~ V M0) = 1. (3.1) 



w p (s — 1) 1 — XIV s ' 

Due to the Fermat curve condition the function w p (s) is cyclic: w p (s + N) = w p (s). 

The explicit formulas for the eigenvectors of open RTC, that is eigenvectors of A n (X), can be 
extracted from [8], where the formulas for the eigenvectors of open Bazhanov-Stroganov quan- 
tum chain (which includes RTC) were obtained. Here we give the resulting formulas specialized 
to open RTC with a little change of notations. 

The eigenvectors of open RTC are given in terms of w p (s) defined by the different points 
Pm,s;m',s' = (x m , s ;m> ,s> , ym,s;m' ,s') , m, vn! = 1, 2, . . . , n; s = 1, 2, . . . , m; s' = 1, 2, . . . , m' , belonging 
to Fermat curve x N + y N = 1. The ^-coordinates of these points are defined by unknown 
parameters (amplitudes) v m ^ s by x mtS - m r >s ' = v m ,sAw,s'- Then the corresponding y m ,s-,m',s' are 
defined up to a root of 1. The parameters v m , s are defined by the initial value v\\ := v and the 
recurrent relations 

Vm\V m 2 ' " ' Vmrn = Vm-l,lV m -l,2 ' ' ' Vm-l,m-lV , 
m 

EI Um-l,l;m,k 

nym—l,s;m—l,l k=l -, 7 -i n -■ 

y — ~ ~—2 = 1 > i = l,2,...,m-l. 

11 !/m-2,s;m-l,i 
s=l 

In the case of homogeneous RTC these relations can be solved explicitly [TT] (see also [8]). 

We will use the notation \y n ) G Vi (8> • • • <£> V n for the eigenvectors of the operator A n (X) of 
the open RTC with n particles. These eigenvectors are labelled by n parameters 7„ jS G Zjy, 
s = 1, 2, . . . , n, collected into a vector 

In = (jn,l, 7n,n) G (Zjv) n . 
The formula 

n n— 1 

l7n) = X] / U (7n-l)Q(7n-l|7n)|7n-l) ® I ^ ~ 

7n-l fe=l Z=l 

where 

n-l 



M(7n-l) = II Wpn-i j!B -i,i(7n-lJ - 7n-l,/ - 1) 



Q(7n-l|7n) 



k/7n-l,lH h7n-l,n-l 



n ra— 1 

n n w p n -i, U n, k (^n-i,i - in,k - 1) 

k=l 1=1 



gives the eigenvectors \-y n ) of .A n (A) from the eigenvectors |7 n -i) £ Vi <8> • • • <8> V n _i of A„_i(A) 
and basis vectors \~/) n G V n . To find the formula for \~fn-i) we can use the same formula and 
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so on. At the last step we need the eigenvectors of 1-particle RTC which are just basis vectors 
|7i,i)i G Vi- The vectors \-y n ) obtained by the described procedure satisfy 



AnWhn) - ft (l + 22^) |7»> = f[ (l - ¥) hw) 

fc=l ^ ' fe=l ^ ' 

B n {\n,k)\ln) = ~^\ln k ), 



n ) ) 



k=ls^k 

n-l 

where = —>tv n ^uJ ln ' k and h k = u \\ y n -l,l;n,k- I n the above formulas the vector |7~ fe ) 

l=i 

means the vector \"f n i, . . . , 7 nn ) with replaced by 7 nj fc — 1. 

We will omit in what follows the index n in matrix elements of monodromy matrix. 

At A = ujfi, the rank of R(X, /i) becomes 1. Therefore both sides of (j2.2[) become proportional 
to i?(A, fj,) with the coefficient being the so-called quantum determinant: 

V 

qdet T n (n) = D{oj^)A(^) - B(u)fi)C(fi) = qdet L^/j.) ■ ■ ■ qdet L n (fi) = — , (3.2) 

{UfJ,) n 

where 

n n 

v = n v k , vi7„) = nK,^ 7n ' fc ) i7»). 
fe=i fc=i 

Now use (|3.2p for \i = \ n y. 

V 



{u\ n ,r) n 

Acting by both sides on |-y n ) we get 



ri 



I] (Vn,k^- k ) 
-B(u\ n>r )C(\ nir )\*1 n ) = — — \j n ). 

Therefore 

n 

(uM nir .) n L h r 

Since C(A) is a polynomial in 1/A of degree n-l we can reconstruct the action of C(A): 

An,r(A An,m) 



— : \ A^A n . r — A nm i 

r=l \m=tr 



Finally using (|3.2p after some calculation we get 



n 

LUX ■ 



X 

sj^k,r 



Vn,k^ n ' k 
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4 Eigenvectors for the one-boundary RTC 

We restrict our attention in what follows to the choice j3± = S± = 0. In this case the boundary 
matrices become 



a-(f-A) r,. 
??_ 



K + {\) 



K-{\) = 
Due to (|2.3p we have 

Z n (A) • • • Li(A) = {uj9) n \~ 1 a z T t {uj9/X)a z = (w#) n V 
Therefore 

A{\) = {ujoyv- 1 



rj + LO 



.1 ( A{uj9/\) -CM/A) 
-B(lo9/X) D(lo9/X) 



9 — A 2 

— : — A(uj9/X)A(X) + r/_(A(w0/A)C(A) - C{wB/\)A[X)) 



X 



B(X) = (cu^V" 1 

c(x) = Mfv- 1 

V(X) = (w0) n V- 1 



x- 



a- 



X 



x- 



-a. 



-a. 



X 

6 - X 2 



X 



A(uj9/X)B(X) + 7]-{A(u6/X)D(X) - C{u9/X)B(X)) 
B(cv9/X)A(X) + 7]^(-B(uj9/X)C(X) + D(u8/\)A(\)) 
B{ujO/\)B{\) + t]_{-B{uj6/X)D(X) + D{uj9 / X)B{X)) 



(4.1) 



From (1221) it follows that A{X) defined by .4(A) = A.4(A)/(0 - A 2 ), or explicitly 



n\r-l 



A(X) = (to9) n V 



a-A(uo6/X)A(X) + 



rj-X 



■{A(uo9/X)C{X) - C{lo9/X)A(X)) 



(4.2) 



(6 - A 2 ) 

satisfies A(<f) = A(X). It means that -4(A) can be presented as 

A{X) = a_x n ((A + Lo9/X) n + (A + uO/X^Hf + • • • + (A + u6/\)H*_i + H*) , (4.3) 

where the set of operators Hf,H2,---,H® have to be identified with the commuting set of 
Hamiltonians for the one-boundary RTC. Explicitly Hf is given by (|5.2p with rj + = 0. Our 
problem is to find the common eigenvectors for this set. 
For a = (cri, . . . , cr n ) G (Z N ) n , let 



^" = ^2Qhn,tT)\ln), 



where the sum includes all N n combinations of -y n = {7n,i ; ■ ■ ■ , 7n,n} and 



r=l fc=lr=l 



(4.4) 



n 



r<r' 



Wp r y{ln,r + Jn,r>) J 
Then, as it will be shown in the Appendix A, we have 

'(A + xs k u ak )(u6/X + xs k bj ak ) 



(4.5) 



^(A)* CT = o_I 



k=l 



S k UJ 
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o_(-x)" II ( (A " Afc)( ^ /A " Afc) ) *«* (4.6) 



k=l 



where Sk, k = 1,2, ... ,n, are some fixed amplitudes which will be defined later. Also we have 
used a short notation Xk = —xskU ak . This form of the eigenvalue polynomial corresponds 
to (|4.3j) . We are using three types of points p r , r ', V^k anc ^ Prk belonging to the Fermat curve 
x N + y N = 1 in the definition of Q(7n, er). They are defined by their ^-coordinates: 

x r y = x 2 v n , r v n y/(u9), x^ k = f re , r /Sfc, xf k = x?v n , r s k / {u)9) 

and the equations 

2 n 

[ y T y = a^h r Y[y® k Vr,k, r = \,2,...,n. (4.7) 



x n-2 a; n+l 



r'^r k=l 



To solve these equations we have to take the iV-th degree of both parts of the equations (|4.7p 
and solve them with respect to elementary symmetric polynomials of s k , k = 1, 2, . . . , n. This 
information is enough to find the values of Sk and the coordinates of the Fermat points entering 

Taking into account Laurent polynomial structure of C(A) in A it is clear that it is enough 
to know the action formulas for C(X) on ^ a in 2n + 2 points of A to reconstruct C(X)^ tT . The 
following formulas are valid: 

C(A fc )M/ (T = r ? _/ fc ^ (T+fc , C{u6/\k)^< T = i 1 -{uO) n f k ~ l ^ (T - k , (4.8) 
C(±Ve)y„ = V _(±Ve) n y„, (4.9) 

n 

where Xk = —xskUJ IJk , fk = Yi ^kljfk/yfk- The reconstructed formula for C(A) is 

r=l 

A2 V-TT Afc(A - \ t )(uj6 - XX r ) 



A ZT k x ( x k - X r )(u9 - A fc A r ) 



+ -YT\ ( A ~ A r)M-AA r ) eX/y/eMey/e)*.. (4.10) 

2 e ^ lr 1 = i 1 A( e ^-A r )(^-A r ) V 1 ^ 

Here we give some heuristic explanation of formulas (14.81) . Let be an eigenvector function 
of A(X) satisfying (|4.6p . Then acting by both sides of commutation relation 

(ufi - n){9 - \p)C{\)A(ji) + {n - X)(uj9 - \fi)A(\)CQi) 

+{fi - X)(l>Xh - \n)B(\)AQi) = {9- Xn)(ufi - X)C(n)A(X), 

which follows from (|2.4p . on ^>a- at n = Xk we get 

Using (|4.6p we see that C(Xk)^a- is an eigenvector of *4(A) with shifted zero, that is Xk replaced 
by ojXk- It means that C(Xk)^er is proportional to ^ a +k. Clearly, this argumentation is not 
sufficient to prove the relations (|4.8p . The derivation of these relations can be carried out 
adapting the corresponding derivation from [5]. To prove (|4.9|) we use (|4.1|) at A = ±y/9 and 
then (D. 
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5 Eigenvectors for the two-boundary RTC. 
Separation of variables 

From (|2.6p it follows 

_ C(X)(co-l)9 + C(f)(6-X 2 ) 

B[X > ~ Ue^x 2 ' 

which gives 

<jj 2 Q — \ 2 

t(X) = a+A(X) + r?+C(A) + r)+u>B(X) 

(^6 - X 2 )(6 - X 2 ) ~ C(A)(^-A 2 )+^C(f)(g-A 2 ) 
= j$ a + A{\) + i] + ^ _ A2 . (5.1) 

We read off the structure of t(A) from this expression 

t(A) = -(u 2 6 2 /X 2 - 9)(X 2 - fl) a+a ~ X " ((A + u>6 / X) n + (A + uj0/X) n ^Hf B + ■■■ 

6 

+ (A + u)0/\)H™ x + O + v+v _ Pn{X ), 

where 

Pn (A) = (A + cu9/X)6 {n - 1)/2 if ?i is odd, p n (X) = (1 + w)# n/2 if n is even 

and the set of operators H® B ,H BB , . . . ,H® B should be identified with the commuting set of 
Hamiltonians for the two-boundary RTC. Explicitly, 

n n—1 
k=l k=l 

- {rj./a^u^vi - (r? + /a + )x~ 1 v~ 1 it n . (5.2) 



Let E = (E\, . . . , E n ), be a common eigenvector for the commuting set of operators 
{^ BB }: 

Hf B $ E = E k $ E 

and <e(A) is corresponding eigenvalue polynomial for t(X). In accordance with a procedure of 
the Separation of Variables we are looking for $e to be in the form 

n 

$E = ]T H((Xi - Xj)(u6 - XiXj)) H q k (a k ,B) 

<t i<j k=l 

where X k = —xs k Lo ak . 

Using the action formula (|5.ip for t(X) on together with ()4.6|) and (j4. 10[) we get Baxter 
equations for separated variables q k (a k ,~E): 



t E (X k )q k (a k ,E) = 77+77, f X *hy - 2n f k q k (<rk ~ 1,E) 

(w0 - w Z A£) 



[V - LOA k ) 
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where fk = Yl ^kyf k /yf k - The corresponding Baxter equation can be considered (for each 

r=l 

k = 1, 2, ... ,n) as a homogeneous system of linear equations for N unknowns q k {a, E), cr G Z^v- 
The requirement of the existence of non-trivial solutions for all k = 1,2, ... , n fixes the common 
spectrum for the commuting set of operators {-ffj? B }. Unfortunately, this problem can not be 
solved explicitly for general N. But if we know the spectrum we can construct the eigenvectors 
solving the systems of linear equations. We would like to report about such a progress [7j in the 
case of the periodic Baxter-Bazhanov-Stroganov model at iV = 2 which contains in particular 
the Ising model. 



6 Another proof of commutativity of 

In this section we describe an alternative approach to the proof of commutativity of Hamiltonians 
using as an example the set of the one-boundary RTC Hamiltonians H?. This approach is based 
on the results from [12]. Let us consider the generation function 

k=l K k=X s^k r s 

= a J] (A + v k + 9/v h ) + r, £ u, 1 V ]J (A + \ + / / ^ ) , (6.1) 

fe=l fe=l s^k Vk s 

where A = A + 6/X, V = v\Vi • • • v n . The set of operators Uj, Vj satisfy the relations 

UjU k = U k Uj, Vji> k = V k Vj, UjV k = q S i- k V k Uj, 

where q is a non-zero complex number. Note that this expression for A(X) coincides with (|4.2ft 
(see also (|A.3P ) up to a redefinition of parameters. 
Expanding A in A we obtain 

A = a (A n + tffA™" 1 + • • • + H^A + H%). (6.2) 

We will prove that the operators obtained by such expansion are pairwise commuting. Let 
us introduce new operators A^ = — v k — 0/v k . Taking into account 

A k - A s = (v k - v s )(6/(v k v s ) - 1) , 

we get 

n n A A 

*-«n<A-A*>+lS>IlCT. 

k=l k=l s^k k s 

where 



n 

Finally using the interpolation formula for the polynomial Y\ (A — A k ) — A" of degree n — 1: 

k=l 

n n A A 

n(A-A t )=A--E«n^ 

k=l k=l s^k K 
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we have 

n A A 



fc=l s ^k " 



where Z^ = r]Zk/a — A r k l . It is easy to verify that 

[A^AfcHO, [Zj,Z k ]=0, [Zj,A fc ]=0, j^k. (6.4) 

Theorem 1. For the set of operators A^, Z\~, k = 1,2, ... ,n, satisfying ()6,4p . we have commu- 
tativity of operators H? defined by expansion ()6.2[) of A given by (|6.3p . 



Proof. Let / be the n by n + 1 matrix with the matrix elements fij, l<i<n,0<j<n, such 
that any two of them belonging to different rows are commuting. Denote by Aj, < j < n, the 
determinant of matrix / with j-th column omitted. Then from [T2] it follows that the operators 
Hk = A/cAq , k = 1, 2, . . . , n, are pairwise commuting. 

Let us fix the matrix elements of / by the matrix elements /^o = Zi, = A™ _fc , where i, 
k = 1, 2, . . . ,n. In this case we have Ao = \\ (A m — A/). Expanding A s , s > 1, with respect to 

m<l 

first column we have 

n 

A s = ]Tz fe (-l) fe+1 A fc; o, s , 
fc=i 

where Afc ; o, s is the determinant of matrix / with fe-th row and 0, s-th columns omitted. Therefore 
for the generating function for commuting H s = A s Aq 1 we get 



it 



fc;0. 



k=l 




The expression in the brackets is the determinant expanded with respect to A:-th row of / with 
0-th column omitted and all A^ replaced by A. Hence this expression is similar to Ao but with A^ 
replaced by A. It gives 

n n A A 

E(-u* +i ".A-«=x>n^. 

s=l k=l s^k K s 

Comparing this formula with (|6.2|) and (|6.3|) we get = {—l) k+1 and therefore the com- 
mutativity of H®. ■ 

Taking into account the discussion before the theorem we get that the function (|6.ip is 
a generating function for the commuting set of operators H®. 



7 Discussion 



We had applied in this contribution the Separation of Variables to obtain the commuting Hamil- 
tonians eigenvectors of the quantum relativistic Toda chain at a root of unity with boundary 
interaction. As we already have discussed in Introduction, the L-operator of the relativistic 
quantum Toda chain at a root of unity is a special case of the cyclic L-operator of the Baxter- 
Bazhanov-Stroganov model (BBS). It was shown in the papers ([3E0])) that the BBS model at 
N = 2 coincides with the free fermion point of the generalized Ising model. We plan to apply 
the Separation of Variables method, used in this contribution, to the N = 2 BBS model with 
integrable boundary conditions. 
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A Derivation of the action formula for -4(A) 

In this Appendix we will prove the formula 

■(\-\ k )(uo/\-\ k y 



A(\)* tr = a-(-x) n l 



k=i 



A* 



(A.l) 



where the eigenvector ^>a- is defined by (14, 4j) . (|4.5p . In what follows we use the short notations 
K,k = -xv n , kW ln ' k and A fc = ->cs k u; ak . 

To prove (|A.ip we use the expression (|4.2p 



-4(A) = {u6) n \- 1 
and the action formulas 



a-A(uB/\)A(\) + 



ry_ A 



(A(lu6/\)C(\) - C(uO/\)A(\)) 



(A.2) 



n / \ \ n 

Mx)hn) = n 1 - ) i7n>, c(A)i 7n ) = 

fc=l ^ ' fc=lr^fc 



An,fc(A — A ni r) 

A(A nj fc — X nr ) 



C{\ n ,k)\ln) 



C{\ n ,k)\ln) 



n 

W A nm 

m=1 |~+*\ 

\n-l I In I 



(->c) n hk(tjj\n,k) r 



We have 



(A{uj9/\)C{\) - <?M/A)4(A))| 7n ) 



A 2 



0\ 



yi *e>* n 

fc=l r^fc 



AA n , 



A / A n — A n . 



C(A n , fc )| 7n >. 



Therefore 



l(A)| 7n > = Mf] «_(-*)« f[ (1-AWM))(1-WA) |7) 



fc=i 



A 



?i. A: 



(A.3) 



fc=l r^fc 



AA n . 



An - 



A / A n t — A n . 



uJ 



hk (wA ni fc) 



n-l 



l7n 



Acting -4(A) on and shifting in an appropriate way the summation we get 
A(X)^ tT = J2Qhn,<T)A(X)\y n ) 



In 



( w\^nt J r v»fr ( 1 - A W M))(l-An, fc /A) 



AA n r 



fc=l 
An,r 



A n ,fc Q(j n k ,cr) 
h k u n+1 Q( 7 n,cr) 



I7n) 



Using (|4T5j) and (|3TTj) and then (14771) we obtain 



Q( 7 - fc ,cr) _ x 2 ^ fc * (A m -A n , fe )(l- A m A n , fc /M)) 



Q( 7 n, O-) A 2 



n 



71, m=l 



AmZ/F m Z/fe 
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rj^k 



U X n ,k — X n , 



X n ,k — <^n,r (1 — ^n,k^n,r / (^0)) 



n,K m=l 



n 

rj^k 



& X n ,k ~ X n ,r 1 



X n ,k — X n>r (1 — ^n,fc^n,r/ {^0)) / 

Finally, taking into account the identity 



n n a ■ n n 

n(A - A„, fc ) + n A r n^™ - a *) = n( A - a *) 

k=l m=l r^m n ' m n ' r k=l k=l 

with 

A = — A — uj9/X, Afc = — Afc — Lo6/Xk, A. n ,k = — X U: k — coO/X n ^ 
we get (TO) . 
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